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Abstract 

A  two-dimensional  barotropic  finite  element  model  with  the  grid  resolution  varying  between  0.2  and  2  km  and  forced 
by  eight  principle  constituents  (Qt,  Ou  Pt,  K\,  N2,  M2 ,  S2,  and  K2)  was  used  to  compute  tidal  elevations  and  currents  in 
the  Bab  el  Mandab  Strait.  Good  agreement  is  achieved  with  the  available  observations  for  both  diurnal  and  semidiurnal 
tidal  currents  and  diurnal  elevations;  however,  the  model  performs  less  satisfactorily  for  the  semidiurnal  elevations 
mainly  due  to  the  errors  between  the  observed  and  computed  phases  in  the  region  where  there  are  amphidromic  points 
for  the  M2,  S2,  and  N2  constituents  in  the  Strait.  The  results  indicate  that  the  largest  amplitudes  of  the  tidal  elevations 
and  the  strongest  currents  are  present  in  the  southern  part  of  the  Strait.  Residual  circulation  induced  by  the  barotropic 
tides  is  rather  weak  in  the  major  part  of  the  Strait,  and  its  contribution  to  the  Red  Sea  water  transport  is  small.  The 
model  results  also  show  that  barotropic  energy  fluxes  are  not  very  large  and  their  direction  depends  on  the  constituent. 
All  diurnal  and  one  semidiurnal  ( N2 )  constituent  have  one  major  source  of  energy,  which  is  the  flux  from  the  Gulf  of 
Aden,  while  there  are  two  sources  of  energy  for  the  M2,  S2,  and  K2  components:  one  from  the  Gulf  and  another  from 
the  Red  Sea.  Very  small  energy  fluxes  from  the  Strait  to  the  adjacent  basins  indicate  that  almost  all  tidal  energy  is 
dissipated  within  the  Strait.  The  distribution  of  the  rate  of  energy  dissipation  due  to  bottom  friction  implies  that  the 
major  area  of  dissipation  is  located  between  Perim  Narrows  and  the  Assab-Mocha  line. 
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1,  Introduction 

The  Bab  el  Mandab  Strait  is  located  in  the 
southern  part  of  the  Red  Sea  and  extends  from 
Perim  Island  in  the  south  to  the  Hanish  Islands  in 
the  northwest  (Fig.  1).  This  Strait  is  the  only 
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Fig.  1.  Map  of  the  Red  Sea,  Bab  el  Mandab  Strait,  Gulf  of  Aden  and  north-western  Indian  Ocean;  shown  are  locations  of  the  water 
level  stations  (stars),  subsurface  pressure  gauges  (triangles),  ADCP  moorings  (closed  circles),  Perim  and  Hanish  transects,  and  open 
ocean  boundary  of  the  model  (dotted  line);  depth  contours  are  in  meters. 


passage  for  the  Red  Sea  waters,  an  important 
component  of  the  thermohaline  circulation  in  the 
Indian  Ocean,  to  the  Gulf  of  Aden/Indian  Ocean 
system. 

A  strong  thermohaline  exchange  flow,  which  is 
the  most  prominent  feature  of  the  circulation 
within  the  Bab  el  Mandab  Strait,  is  highly 
influenced  by  both  synoptic  timescale  wind  forcing 
and  strong  tidal  motions.  This  exchange  flow 
exhibits  a  distinct  seasonal  variability,  and  the 
overall  picture  emerging  from  all  available  inves¬ 
tigations  is  a  two-layer  winter  inverse  estuarine 
(fresh  inflow  on  top  of  a  hypersaline  outflow) 
exchange  during  October  to  May,  replaced  be¬ 
tween  June  and  September  by  a  three-layer 
exchange  comprised  of  a  shallow  surface  outflow. 


an  intermediate  intrusion  of  the  relatively  fresh 
and  cold  Gulf  of  Aden  Intermediate  Waters,  and  a 
deep  hypersaline  outflow  which  is  a  small  fraction 
of  the  winter  value  (Vercelli,  1931;  Thompson, 
1939;  Smeed,  1997;  Murray  and  Johns,  1997). 
Recent  and  historical  observations  indicate  that 
tidal  motion  is  also  a  very  significant  factor  of  the 
flow  field  in  the  Bab  el  Mandab  Strait  (Vercelli, 
1925,  1927;  Siedler,  1969;  Jarosz,  2002).  The  tidal 
elevations  and  currents  there  exhibit  strong 
diurnal  and  semidiurnal  oscillations.  Barotropic 
tides  dominate  in  this  Strait,  while  baroclinic  tides, 
which  are  more  prominent  for  the  winter  exchange 
flow  period  and  primarily  of  a  diurnal  period,  are 
less  energetic  (Jarosz,  2002).  The  strongest  internal 
tidal  fluctuations  are  confined  to  a  pycnocline 
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layer,  i.e.  to  depths  within  the  water  column  where 
the  greatest  density  gradients  are  observed.  In 
addition,  analyses  of  historical  data  led  to  the 
conclusion  that  there  is  an  amphidrome  for  the  M2 
tidal  wave  in  the  Bab  el  Mandab  Strait  (Vercelli, 
1925;  Defant,  1961). 

The  available  observations,  however,  are  spa¬ 
tially  scattered  and  too  sparse  to  examine  tidal 
dynamics  in  the  Bab  el  Mandab  Strait  in  more 
detail.  It  is  rather  difficult  with  these  data  to 
conclude,  for  instance,  whether  there  are  also 
amphidromic  points  for  other  major  tidal  consti¬ 
tuents  and  where  these  points  are  located.  With 
this  data  set,  it  is  also  difficult  to  evaluate  tidal 
residual  currents,  their  variability  and  possible 
contribution  to  the  transport  of  the  Red  Sea 
waters  to  the  Gulf  of  Aden/Indian  Ocean  system 
or  to  estimate  barotropic  tidal  energy.  To  better 
understand  the  tidal  dynamics  in  the  Bab  el 
Mandab  Strait  and  address  some  of  the  posed 
questions,  a  numerical  model  with  a  very  high 
spatial  resolution  is  implemented.  Specific  objec¬ 
tives  are  to  (1)  quantify  the  spatial  variability  of 
the  barotropic  tidal  elevations  and  currents  in  the 
Bab  el  Mandab  Strait,  (2)  evaluate  the  barotropic 
tidal  residual  circulation  and  its  contribution  to 
the  transport  of  the  Red  Sea  waters,  and  (3) 
estimate  the  barotropic  energy  fluxes  and  dissipa¬ 
tion  rates  for  this  Strait. 

The  barotropic  tides  in  the  Strait  are  numeri¬ 
cally  simulated  with  a  two-dimensional  form 
of  the  finite  element  hydrodynamic  model 
called  ADvanced  CIRCulation  Model  for  Shelves, 
Coasts  and  Estuaries  (ADCIRC)  developed  by 
Luettich  et  al.  (1992)  and  Westerink  et  al.  (1994). 
Selection  of  a  finite  element  model  is  justified 
by  its  grid  flexibility,  which  allows  easy  incorpora¬ 
tion  of  the  complex  coastline  and  high  levels  of 
refinement  near  shallow  coastal  areas  and  in 
regions  of  rapid  bathymetric  change.  It  also 
permits  coarsening  of  the  grid  over  deeper 
waters  to  extend  the  domain  boundaries 
away  from  the  primary  region  of  interest. 
Furthermore,  finite  element  models  have  a  history 
of  accuracy  in  modeling  barotropic  tides  (Walters, 
1987;  Werner  and  Lynch,  1987;  Walters  and 
Werner,  1989;  Westerink  et  al.,  1989,  1992;  Fore¬ 
man  et  al.,  1995). 


This  paper  is  organized  as  follows.  In  Section  2, 
the  numerical  model  and  finite  element  mesh  is 
briefly  described.  Model-data  comparisons  are 
presented  in  Section  3.  Variability  of  tidal  eleva¬ 
tions  and  currents  are  discussed  in  Section  4. 
Residual  circulation  is  described  in  Section  5,  with 
energy  fluxes  and  dissipation  rate  distribution 
presented  in  Section  6.  Section  7  summarizes  the 
findings  of  the  study. 


2.  Model  description 

The  model  domain,  shown  in  Fig.  1,  includes 
not  only  the  Bab  el  Mandab  Strait  but  also  the 
entire  Red  Sea,  Gulf  of  Aden  and  northwestern 
part  of  the  Indian  Ocean.  Such  a  large  model  area 
was  chosen  primarily  to  reproduce  tidal  waves 
propagating  from  the  Indian  Ocean,  which  is  a 
major  forcing  of  tidal  motion  in  the  Strait  and  Red 
Sea  as  discussed  by  Defant  (1961).  A  second 
motivation  was  to  avoid  a  necessity  of  having  two 
open  boundaries  at  the  northern  and  southern 
ends  of  the  Strait.  The  tidal  elevation  data  outside 
and  within  the  Strait  are  limited,  thus  estimation 
of  the  reliable  tidal  elevation  harmonic  constants 
at  the  open  boundaries  located  near  the  Strait 
would  be  difficult. 

Bathymetry  for  the  model  was  obtained  from 
two  sources:  the  Naval  Oceanographic  Office 
Digital  Bathymetric  Data  Base — Variable  Resolu¬ 
tion  (DBDB-V)  (NAVOCEANO,  1997)  and  charts 
published  by  the  Defense  Mapping  Agency  in 
1992.  The  bathymetry  from  the  charts  was 
digitized,  and  this  data  set  (of  resolution  ~0.11') 
is  limited  to  an  area  slightly  larger  than  the  Strait 
itself,  while  the  DBDB-V  data  with  a  resolution  of 
2'  cover  the  remainder  of  the  model  domain. 

The  finite  element  grid  used  in  computations  is 
displayed  in  Fig.  2.  It  consists  of  37,436  nodes  and 
70,733  elements.  Nodal  spacing  for  this  mesh 
varies  throughout  the  modeled  region  and  ranges 
between  0.2  and  55.5  km  with  the  highest  refine¬ 
ment  present  in  the  Strait  where  the  minimum  and 
maximum  nodal  spacing  are  0.2  and  2  km, 
respectively.  The  coarsest  resolution  is  located  in 
deep  waters  of  the  Gulf  of  Aden  and  Indian 
Ocean. 
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Fig.  2.  Model  finite  element  grid. 


The  two-dimensional  form  of  the  ADCIRC  model 
is  based  on  vertically  integrated  equations  of  motion 
and  continuity,  which,  in  a  spherical  coordinate 
system,  are  defined  as  follows  (Gill,  1982): 


S  U  U  0  U  VdU  UV  sin  tp 

3 1  Rcostp  02  R  dtp  Rcostp 

9  d  Tbx 

RcostpdX^  at])  p0H 


(la) 


dv  U  dV  V  dV  UV  sin  0 

dt  R  cos  tp  02  R  dtp  ^  R  cos  tp 


9  6  ,r  \  ^b(p 


(lb) 


is  the  total  water  column  depth,  h  is  the  bathymetric 
depth  relative  to  the  geoid,  f  —  2Q  sin  tp  is  the 
Coriolis  parameter,  Q  is  the  angular  speed  of  the 
Earth,  p0  is  a  reference  density,  g  is  the  acceleration 
due  to  gravity,  a  is  the  Earth  elasticity  factor 
approximated  as  0.69  for  all  tidal  constituents  as 
used  by  other  investigators  including  Schwiderski 
(1980)  and  Hendershott  (1981);  however,  its  value 
has  been  shown  to  be  constituent-dependent  (Wahr, 
1981),  tj  is  the  Newtonian  equilibrium  tidal  potential, 
and  r bx,  Tbcp  are  the  bottom  stresses  taken  as: 


m  =  P0cduVu2  +  v2, 
xb<t,  =  p0cAvVu2+v2, 


0C  1  8 (UH)  1  d(  VH  cos  tp)  Q 

dt^  Rcostp  02  ^  Rcostp  dtp 

(lc) 


where  Cd  denotes  the  bottom  drag  coefficient.  The 
equilibrium  tidal  potential  is  expressed  as  (Reid, 
1990) 


where  t  represents  time,  A,  tp  denote  degrees  of 
longitude  and  latitude,  '(  is  the  free  surface  elevation, 
U,  V  are  the  depth-averaged  horizontal  east-west 
and  north-south  velocities,  respectively,  H  =  £  +  h 


*7(2,  (p,t)  =  Y1  cjnf  jn(<b)Lj(fp)  cos  [271  (t  -  t0)/Tjn 
nj 

+jX  +  ^(/o)], 


(3) 
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where  t  is  time  relative  to  to,  which  is  the  reference 
time,  Cjn  is  a  constant  characterizing  the  amplitude 
of  a  tidal  constituent  n  of  species  j,  fjn  is  the  time- 
dependent  nodal  factor,  ojn  is  the  time-dependent 
astronomical  argument,  j  =  0,  1,  2  are  the  tidal 
species  (j  =  0  declinational;  j  =  1  diurnal,  j  =  2 
semidiurnal),  Lo  =  3  sin2  cp  -  1,  L\  —  sin  (2 cp),  L2  = 
cos2  cp ,  and  Tjn  is  the  period  of  a  constituent  n  for 
species  j.  The  tidal  loading  and  self-attraction  are  not 
included  in  the  model  equations. 

Prior  to  being  discretized,  the  continuity  and 
momentum  equations  are  combined  into  a  general¬ 
ized  wave  continuity  equation  (GWCE),  which  has 
been  shown  to  have  superior  numerical  properties 
to  a  primitive  continuity  equation  when  a  finite 
element  method  is  used  in  space  (Lynch  and  Gray, 
1979).  The  final  forms  of  the  GWCE  and 
momentum  equations,  which  are  solved  by  the 
model,  are  given  in  Blain  and  Rogers  (1998);  while 
numerical  discretization  of  these  equations  is 
described  in  detail  by  Luettich  et  al.  (1992)  and 
Kolar  et  al.  (1994). 

No-normal  flow  and  free  tangential  slip  condi¬ 
tions  are  imposed  at  the  land  boundary  nodes. 
Furthermore,  at  the  open  ocean  boundary,  the 
tidal  elevation  generated  by  four  diurnal  (A),  0\, 
P i,  Qi)  and  four  semidiurnal  ( M2 ,  S2,  N2,  K2) 
constituents  is  specified  (Dirichlet  condition).  The 
tidal  harmonic  constants  used  to  generate  the 
elevation  along  the  open  boundary  are  linearly 
interpolated  onto  the  boundary  nodes  using  data 
from  the  World  Ocean  Tide  Model  database 
FES95.2  (Le  Provost  et  al.,  1994).  In  addition, 
an  equilibrium  tidal  potential  forcing  within  the 
domain  is  applied  for  the  same  eight  constituents. 

A  constant  value  for  the  bottom  friction 
coefficient,  equal  to  0.003,  was  applied  throughout 
the  domain.  Several  bottom  friction  coefficients 
were  tried,  and  the  value  of  0.003  gave  the  best 
agreement  with  the  observations.  A  time  step  of 
30  s  was  used  to  ensure  model  stability  based  on 
the  Courant  number  criterion.  The  parameter  to, 
which  weights  the  primitive  and  GWCE  form  of 
the  continuity  equation,  was  estimated  from  a 
formula  given  by  Westerink  et  al.  (1994)  and  set 
equal  to  0.001.  Finally,  the  minimum  depth  was 
assigned  to  be  2  m  to  eliminate  any  potential 
drying  of  computational  nodes  since  wetting  and 


drying  option  is  not  included  in  these  simulations. 
The  model  simulations  were  carried  out  for  1  year 
to  generate  a  long  time  series  that  allows  the 
separation  of  P\  and  K2  constituents  from  K\  and 
S2 ,  respectively.  The  amplitudes  and  phases  of  the 
tidal  constituents  were  obtained  through  the 
standard  harmonic  analysis  (Foreman,  1977, 
1978). 


3.  Validation  of  computed  tides 

Tidal  harmonic  constants  at  29  water  level 
stations  and  four  subsurface  pressure  gauges  as 
well  as  estimates  of  the  barotropic  tidal  currents 
obtained  from  the  measurements  at  five  ADCP 
moorings  collected  during  the  project  entitled 
“Observation  and  Modeling — an  Integrated  Study 
of  the  Transport  through  the  Strait  of  Bab  el 
Mandab”  (the  BAM  project)  (Murray  and  Johns, 
1997)  were  used  to  evaluate  computed  tidal 
elevations  and  currents.  The  ADCP  observations 
were  vertically  averaged  before  they  were  com¬ 
pared  to  the  modeled  currents  (see  Jarosz  (2002) 
for  more  details  regarding  pressure  gauge  and 
ADCP  data  processing  and  analyses).  Geographi¬ 
cal  locations  of  the  water  level  stations,  for  which 
only  the  tidal  harmonic  constants  were  obtained 
from  the  International  Hydrographic  Organiza¬ 
tion  (1979),  and  pressure  gauges,  deployed  for  the 
BAM  project,  are  listed  in  Table  1  and  displayed  in 
Fig.  1.  All  observed  tidal  amplitudes  and  phases 
come  from  near-coastal  areas  leaving  no  possibi¬ 
lity  of  verifying  the  model  predictions  of  the  tidal 
elevations  in  open  waters.  In  addition,  at  some 
water  level  stations  only  a  few  tidal  constituents 
were  available  for  the  comparison;  the  common 
constituents  for  all  stations  were  K\,  0\,  M2,  and 
S2. 

Major  features  of  the  diurnal  tidal  elevations, 
such  as  larger  amplitudes  in  the  Gulf  of  Aden  and 
much  smaller  ones  in  the  Red  Sea  proper,  are 
qualitatively  captured  by  the  model  when  com¬ 
pared  to  the  available  observations  (Table  1). 
Similarly,  well  represented  are  major  features  of 
the  semidiurnal  tidal  elevations.  Such  semidiurnal 
features  include  higher  amplitudes  in  the  Gulf 
of  Aden,  Aquaba  and  Suez,  very  small  M2 
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Table  1 

Comparison  of  computed  (Com)  and  observed  (Obs)  AT,/0,  and  M2/S2  amplitudes  (A)  and  phases  (g,  GMT)  at  the  coastal  stations 
Station  name  Latitude(N)/Longitude(E)  ATj/Oj  M2!S2 

A  (cm)  g  O  A  (cm)  g  (°) 

Obs  Com  Obs  Com  Obs  Com  Obs  Com 


Aden” 

12°47' 

40.0 

39.0 

350 

346 

48.0 

49.5 

134 

133 

44°59' 

20.0 

20.0 

352 

350 

21.0 

21.8 

159 

148 

Aqabab 

29°31' 

2.0 

1.7 

158 

160 

28.0 

16.2 

128 

134 

35°00' 

1.0 

0.5 

146 

164 

8.0 

6.2 

155 

154 

Ashrafi  Islands’3 

27o4r 

2.0 

2.0 

167 

159 

25.0 

14.2 

118 

129 

33°43' 

1.0 

0.7 

153 

154 

4.0 

5.1 

145 

151 

Assabc 

13°00' 

18.0 

14.5 

335 

340 

6.9 

14.3 

259 

118 

42°44' 

8.5 

6.7 

344 

344 

4.0 

7.1 

170 

155 

Berbera" 

10°26' 

46.0 

38.3 

349 

348 

48.0 

49.5 

135 

134 

45°00' 

19.0 

19.6 

356 

351 

20.0 

21.8 

160 

161 

Djibouti" 

1 1°35' 

39.0 

38.7 

354 

347 

46.3 

51.4 

139 

134 

43°09' 

19.0 

19.8 

357 

351 

20.5 

22.0 

163 

160 

G14a 

12°38.5' 

39.0 

38.4 

340 

348 

47.0 

50.2 

119 

133 

43°54.5' 

20.0 

19.7 

344 

351 

20.0 

22.0 

139 

160 

G89c 

12°43.5' 

30.0 

27.9 

340 

346 

23.0 

34.4 

121 

131 

43°08' 

15.0 

14.2 

345 

349 

14.0 

15.6 

142 

160 

G109c 

12°43.6' 

30.0 

30.0 

340 

350 

29.0 

34.0 

125 

134 

43°28' 

15.0 

15.2 

345 

354 

16.0 

16.7 

144 

162 

G108c 

13°40.5' 

6.0 

4.4 

321 

350 

24.0 

8.0 

286 

350 

42°  10. 5' 

2.0 

1.5 

335 

354 

5.0 

2.0 

299 

15 

Harmil  Islandb 

16°29' 

2.0 

2.0 

166 

161 

13.0 

15.0 

318 

318 

40°1 1' 

1.0 

0.7 

180 

160 

3.0 

6.4 

334 

344 

Hudaidab 

14°50' 

1.0 

1.2 

340 

27 

30.0 

13.0 

305 

320 

42°56' 

1.0 

0.6 

92 

104 

6.0 

4.5 

351 

341 

Jeddahb 

21°31' 

2.8 

3.0 

156 

159 

6.0 

3.4 

109 

149 

39°08' 

1.0 

2,0 

161 

159 

1.0 

1.4 

132 

172 

Kamaranb 

1 5°20' 

2.0 

1.0 

34 

99 

33.0 

18.0 

300 

324 

42°36' 

1.0 

1.0 

140 

136 

9.0 

6.9 

334 

348 

Marabatb 

16°59' 

33.0 

34.5 

345 

342 

32.0 

32.3 

144 

144 

54°41' 

19.0 

18.2 

346 

345 

14.0 

12.9 

170 

169 

Massawab 

15°37' 

2.3 

2.8 

164 

166 

33.4 

34.4 

328 

333 

39°28' 

2.0 

2.1 

184 

163 

12.4 

15.0 

332 

3 

Mochac 

I3019, 

7.0 

16.0 

335 

352 

8.0 

14.0 

244 

134 

43-14' 

6.1 

7.8 

352 

358 

4.5 

7.2 

188 

166 

Muhammadb 

20°54' 

3.0 

2.9 

166 

160 

6.0 

2.9 

132 

176 

37°  10' 

2.0 

1.7 

175 

159 

1.0 

1.2 

185 

210 

Mulalla“ 

14°32' 

40.0 

36.8 

349 

345 

40.0 

38.6 

136 

134 

49°08' 

20.0 

19.2 

352 

348 

12.0 

16.6 

164 

159 

Perimc 

12°38' 

35.0 

34.5 

350 

350 

37.0 

42.6 

136 

136 

43°24' 

18.0 

18.1 

351 

353 

17.0 

18.9 

159 

162 

Port  Salalah" 

16°56' 

36.0 

34.7 

344 

343 

31.0 

31.9 

144 

143 

54°00' 

18.0 

18.2 

347 

345 

12.0 

12.8 

168 

168 

Port  Sudanb 

19°36' 

2.0 

2.9 

168 

157 

1.0 

1.5 

204 

246 

37°  14' 

2.0 

1.7 

170 

156 

1.0 

0.8 

256 

277 

Quseirb 

26°06' 

2.0 

2.2 

158 

154 

22.0 

13.0 

112 

134 

34°  16' 

2.0 

1.0 

192 

153 

5.0 

4.9 

139 

156 

Ras  Alula“ 

1 1°59' 

37.1 

35.6 

348 

356 

36.3 

36.5 

139 

138 

50°47' 

18.9 

18.5 

353 

350 

15.4 

15.8 

156 

163 

Ras  Ghan’db 

28°21' 

2.0 

2.9 

160 

163 

18.0 

9.0 

274 

302 
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Table  1  ( continued ) 


Station  name 

Latitude(N)/Longitude(E) 

K\/0] 

m2/s2 

A  (cm) 

9  0 

A  (cm) 

9  0 

Obs 

Com 

Obs 

Com 

Obs 

Com 

Obs 

Com 

33°07' 

2.0 

0.9 

157 

154 

7.0 

3.0 

302 

318 

Ras  Khathibb 

14°55' 

4.0 

1.2 

69 

25 

26.0 

13.0 

294 

318 

42°54' 

1.0 

0.6 

82 

82 

7.0 

4.5 

339 

339 

Saylac“ 

11°22' 

40.0 

38.9 

348 

348 

50.0 

51.6 

137 

133 

43°28' 

20.0 

19.9 

356 

351 

20.0 

22.6 

153 

161 

Shaker  Islandb 

27°27' 

2.0 

2.0 

167 

154 

25.0 

14.0 

117 

132 

34°02' 

1.0 

0.7 

178 

151 

4.0 

5.1 

144 

154 

Sherm  Rabeghb 

22°44' 

4.0 

2.7 

156 

160 

11.0 

6.1 

124 

142 

38°58' 

4.0 

1.4 

162 

163 

2.0 

2.5 

165 

164 

Siqa 

12°4( Y 

35.0 

33.0 

338 

346 

23.0 

27.0 

139 

140 

54°04' 

17.0 

17.3 

339 

349 

13.0 

11.4 

156 

164 

Suez*5 

29°56' 

4.5 

3.3 

158 

162 

56.0 

26.0 

278 

300 

32°33' 

1.3 

0.8 

170 

152 

14.0 

8.7 

306 

319 

Torb 

28°  14' 

4.0 

2.6 

164 

157 

8.0 

3.2 

205 

213 

33037, 

2.0 

0.8 

159 

152 

1.0 

1.6 

230 

211 

Zafaranab 

29°07' 

3.0 

3.3 

165 

163 

42.0 

21.0 

280 

300 

32°40' 

1.0 

0.9 

199 

152 

12.7 

6.9 

301 

320 

“Stations  located  in  the  Gulf  of  Aden  and  Indian  Ocean. 
bStations  located  in  the  Red  Sea. 

“Stations  located  in  the  Bab  el  Mandab  Strait. 


amplitudes  near  Port  Sudan  and  Jeddah  resulting 
from  the  presence  of  the  anticlockwise  amphidro- 
mic  system,  and  the  M2  amphidrome  in  the  Bab  el 
Mandab  Strait  (Defant,  1961). 

A  measure  of  the  agreement  between  the 
modeled  and  observed  elevations  is  achieved  by 
examining  correlation  coefficients  between  com¬ 
puted  and  measured  amplitudes  and  phases.  Table 
2  lists  these  coefficients  together  with  the  number 
of  stations  used  in  the  comparison  and  the 
respective  standard  deviations  of  the  differences 
between  the  modeled  and  observed  values.  In 
addition,  the  parameter  H is  presented,  which  is  an 
average  difference  between  the  observations  and 
model  solution  and  measures  overall  model 
performance.  This  parameter  is  estimated  from 
the  following  expression  (Davies  et  al.,  1997): 

1  N 

H  =  TC  ^  ^[(^obs  COS  t70(-,s  —  ■'^com  COS  gcom ) 
iv  k=] 

T  (•''lobs  sin  gobs  ^com  sinfifcom)2]1/2,  (4) 


Table  2 

Correlation  coefficients  (Corr)  and  standard  deviations  (SD) 
between  the  computed  and  measured  tidal  amplitudes  and 
phases 

Tidal  constituent  N  Amplitude  Phase  H  (cm) 

Corr  SD(cm)  Corr  SD(deg) 


Q\ 

10 

0.97 

0.4 

0.90 

11 

0.53 

0 1 

33 

0.99 

0.8 

0.99 

16 

1.15 

A 

17 

0.99 

0.6 

0.98 

22 

0.93 

K\ 

33 

0.99 

2.4 

0.86 

17 

2.45 

n2 

23 

0.73 

3.2 

0.85 

42 

4.49 

m2 

33 

0.83 

9.3 

0.85 

39 

9.39 

S2 

33 

0.94 

2.5 

0.79 

18 

2.99 

K2 

20 

0.91 

1.1 

0.75 

28 

1.68 

N  is  the  number  of  the  stations;  H  is  an  average  difference 
between  the  observations  and  computed  solution. 

where  N  is  the  number  of  the  water  level  stations,  A 
and  g  are  amplitudes  and  phases,  respectively,  and 
suffixes  “com”  and  “obs”  denote  the  computed  and 
observed  harmonic  constants,  respectively. 
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High  correlation  coefficients  (greater  than  0.96 
for  the  amplitudes  and  0.85  or  higher  for  the 
phases)  coupled  with  the  small  H  values  (less  than 
2.5  cm)  and  low  standard  deviations  (less  than 
2.5  cm  for  amplitudes  and  23°  for  phases)  indicate 
accurate  replication  of  the  diurnal  components  by 
the  model.  For  the  major  semidiurnal  constituents, 
the  model  computations  do  not  compare  as 
favorably  to  the  observations  as  shown  by  the 
lower  correlation  coefficients  (less  than  0.95  for 
amplitudes  and  below  0.85  for  phases)  and  higher 
values  of  H  (~3  cm  or  more  for  the  S2,  N2,  and  M2 
components).  Note  (Table  1  and  Fig.  3)  that  the 
larger  errors  of  the  semidiurnal  elevation  are  not 
uniformly  distributed  within  the  model  domain. 


Fig.  3  shows  mean  absolute  amplitude  and  phase 
differences  computed  separately  for  the  stations 
located  in  the  Red  Sea  (17  stations),  Bab  el 
Mandab  Strait  (six  stations)  and  Gulf  of  Aden/ 
Indian  Ocean  (10  stations).  Comparison  of  the 
measured  data  to  the  model  solution  (Table  1) 
indicate  that  in  the  Gulf  of  Aden  and  southern 
part  of  the  Bab  el  Mandab  Strait,  the  agreement 
between  the  observations  and  the  model  is  very 
satisfactory  for  both  diurnal  and  semidiurnal 
components.  The  same  conclusion  can  be  reached 
for  the  Gulf  Aden/Indian  Ocean  stations  displayed 
in  Fig.  3  by  examining  the  means,  which  are  small 
for  both  amplitude  (less  than  2  cm)  and  phase  (less 
than  7°).  This  good  agreement  suggests  that  the 


Red  Sea  Bab  e!  Mandab  G.of  Aden/Indian  Q 


(d) 


’  95%  Confidence  Interval  for  Mean  Amplitude  Difference] 
Mean  Amplitude  Difference 

95%  Confidence  Interval  for  Mean  Phase  Difference 
Mean  Phase  Difference 


1 7  i  i 


i 


H  125 


100 


150 


75  E 


50 


25 


Red  Sea 


Bab  el  Mandab  G.of  Aden/Indian  O. 


Fig.  3.  Mean  absolute  amplitude  and  phase  differences  between  the  observations  and  modeled  solution  and  their  95%  confidence 
intervals. 
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tidal  elevation  constants  taken  from  the  World 
Ocean  Tide  Model  (Le  Provost  et  al.,  1994) 
represent  well  the  tidal  wave  incoming  from  the 
Indian  Ocean.  Therefore,  the  open  ocean  bound¬ 
ary  forcing  is  not  likely  to  be  a  major  source  of  the 
errors  observed  north  of  Perim  Narrows  where  the 
semidiurnal  waves  are  not  well  simulated,  espe¬ 
cially  with  regard  to  their  phases.  The  largest 
phase  discrepancies  are  found  in  the  Strait  near  the 
amphidrome  for  the  M2  (Assab,  G108  and  Mocha) 
where  the  difference  between  the  computed  and 
observed  phases  can  be  as  large  as  140°. 

A  set  of  experiments,  in  which  alternate  bathy¬ 
metric  values  were  assigned  at  the  nodes  within  the 
Strait,  indicates  that  the  phase  errors  are  partly 
related  to  the  bathymetry.  A  10  m  decrease  of 
depth  in  the  deep  channel  and  a  5m  decrease  of 
depth  in  the  shoals  (at  mesh  nodes  with  depths  of 
7  m  or  more)  between  Perim  Narrows  and  Assab- 
Mocha  line  resulted  in  a  15%  increase  of  the  M2 
phase,  reducing  the  phase  error  and  causing  small 
changes  in  the  K\  phase.  An  increase  of  the  bottom 
friction  coefficient  (Cf)  to  values  of  0.05  or  larger 
for  the  original  bathymetric  configuration  had 
similar  effects  on  the  M2  phase  error,  i.e.  a 
reduction.  Additionally,  it  was  shown  by  Grenier 
et  al.  (1995)  that  in  regions  with  strong  rectilinear 
currents,  the  standard  quadratic  formulation  for 
the  friction  term  might  not  be  sufficient  to 
properly  describe  damping  of  the  tidal  waves. 
Since  the  tidal  currents  in  the  Strait  are  quite 
strong  and  very  rectilinear,  the  quadratic  bottom 
friction  formulation  could  be  introducing  an 
additional  error  into  the  computations.  This  error, 
however,  is  likely  small  in  comparison  to  errors 
associated  with  other  forms  of  tidal  dissipation 
such  as  form  drag,  turbulence,  internal  tide 
generation,  and  interaction  with  nontidal  currents, 
which  are  not  represented  in  the  two-dimensional 
barotropic  model. 

Locations  (see  Fig.  1  and  Table  3  for  the 
mooring  locations)  with  observational  estimates 
(5)  of  the  barotropic  tidal  currents  are  limited  to 
the  Bab  el  Mandab  Strait.  In  Table  3,  the  observed 
and  calculated  Ku  0\,  M2,  and  S2  tidal  currents 
are  described  by  the  parameters  of  the  tidal  ellipse 
(Pugh,  1987).  In  general,  the  agreement  between 
computed  and  measured  tidal  currents  in  the  Strait 


is  good  for  both  semidiurnal  and  diurnal  consti¬ 
tuents.  The  largest  error  for  the  semimajor  axis  is 
found  for  the  M2  at  Perim  Narrows  (moorings 
A2b,  AA1,  and  AA2b)  where  the  computed 
amplitudes  are,  on  average,  20%  smaller  than 
the  measured  ones.  For  other  constituents  at 
Perim  Narrows  as  well  as  for  all  tidal  components 
at  C  mooring  and  the  Hanish  Sill  (mooring  B2b) 
the  difference  between  computed  and  observed 
semimajor  axes  is  always  less  than  20%  of  those 
estimated  solely  from  the  observations.  The 
modeled  semiminor  axes  are  very  small  (2  cm/s 
or  less)  and  comparable  with  those  computed  from 
the  observations.  The  difference  between  the 
model  and  observed  inclinations  is,  on  average, 
4%  of  the  observed  inclination  angles  (measured 
counterclockwise  from  east),  while  the  phase 
difference  is  always  smaller  than  1  h  for  all 
considered  constituents.  Additionally,  the  H  va¬ 
lues  computed  separately  for  along  and  cross-strait 
velocity  components  are  listed  in  Table  4.  The  H 
values,  except  the  M2  component,  are  small 
(<2.30 cm/s),  indicating  that  in  the  Strait,  the 
barotropic  tidal  currents  are  well  simulated  by  the 
model. 


4.  Tidal  elevations  and  currents  in  the  Strait 

Fig.  4  displays  computed  coamplitudes  (in  cm) 
and  cophases  (in  degrees,  GMT)  of  the  tidal 
elevations  for  the  Ku  0\,  M2,  S2,  N2,  and  K2 
constituents.  Similar  to  findings  from  data  ana¬ 
lyses  (Vercelli,  1925,  1927;  Siedler,  1969;  Jarosz, 
2002),  K\  (0j,  Pj  and  Qx  amplitudes  are,  on 
average,  50%,  31%,  and  9%,  respectively,  of  the 
K\)  and  M2  (amplitudes  of  the  S2,  N2  and  K2  are, 
on  average,  45%,  48%,  and  9%,  respectively,  of 
the  M2 )  are  the  dominant  constituents  in  the 
Strait.  Furthermore,  the  amplitudes  of  all  tidal 
components  show  very  similar  behavior  in  the 
Strait,  i.e.  little  variability  in  the  cross-strait 
direction  and  rapid  attenuation  of  the  amplitudes 
in  the  along-strait  direction.  The  largest  values  are 
found  near  Perim  Narrows  and  the  smallest  ones 
are  present  at  the  northern  end  of  the  Strait.  The 
amplitude  and  phase  distributions  indicate  that  the 
M2  and  S2  waves  have  amphidromic  systems 
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Table  3 

Comparison  between  observed  (Obs)  and  computed  (Com)  tidal  current  ellipse  parameters:  semimajor  axis  ( M ),  semiminor  axis  (Mn), 
inclination  angle  (0,  measured  counterclockwise  from  east),  and  phase  ( g ,  GMT) 


Mooring 

Latitude(N)/Longitude(E) 

Ellipse 

Parameters 

0 , 

m2 

52 

Obs 

Com 

Obs 

Com 

Obs 

Com 

Obs 

Com 

A2b 

12°45.12' 

M  (cm/s) 

29.7 

29.4 

16.4 

17.3 

29.2 

25.0 

10.2 

9.1 

43°16.89' 

Mn(cm/s)a 

-1.5 

-1.4 

-0.5 

-1.0 

-1.0 

-0.7 

-0.3 

-0.2 

6  O 

91.3 

98.8 

92.5 

97.9 

94.3 

99.0 

94.3 

99.1 

an 

68.3 

73.2 

68.5 

74.7 

200.8 

221.0 

227.0 

246.0 

AA1 

12°43.26' 

M  (cm/s) 

31.9 

30.0 

18.7 

17.4 

32.5 

24.6 

11.2 

9.5 

43°14.61' 

Mn(cm/s)a 

-0.9 

-1.7 

-0.1 

-1.0 

-0.9 

-0.7 

-0.6 

-0.2 

6  0 

104.4 

104.6 

104.3 

104.6 

105.4 

104.6 

105.2 

104.3 

a  n 

88.6 

72.1 

89.2 

75.0 

234.0 

222.3 

258.5 

246.5 

AA2b 

12°43.64' 

M  (cm/s) 

29.0 

28.9 

16.7 

16.1 

29.6 

23.9 

9.3 

8.8 

43°  16.60' 

Mn(cm/s)a 

-0.7 

-2.0 

-0.3 

-0.9 

-0.5 

-0.4 

-0.2 

-0.1 

6  0 

94.6 

95.0 

94.2 

94.3 

96.9 

96.6 

96.9 

96.0 

9  0 

71.3 

72.1 

74.6 

73.9 

200.8 

220.0 

223.4 

246.0 

C 

13°20.48' 

M  (cm/s) 

13.8 

15.1 

7.2 

8.6 

16.2 

13.4 

6.2 

5.5 

42°54.22' 

Mn(cm/s)a 

1.2 

0.6 

0.6 

0.4 

0.2 

0.3 

0.2 

0.1 

6  0 

133.5 

132.8 

132.0 

133.3 

128.3 

132.7 

124.7 

132.0 

aO 

73.4 

74.6 

74.4 

76.6 

220.9 

222.9 

248.1 

248.9 

B2b 

13°42.62' 

M  (cm/s) 

15.6 

16.5 

9.4 

9.1 

15.0 

15.0 

6.0 

6.4 

42°32.35' 

Mn(cm/s)° 

-0.6 

-0.4 

-0.4 

-0.3 

-0.3 

-0.3 

-0.1 

-0.1 

6  O 

133.0 

122.0 

134.2 

123.0 

133.9 

123.0 

133.4 

123.8 

9  O 

73.9 

78.3 

77.6 

81.8 

204.6 

223.0 

230.3 

250.0 

aA  sign  of  the  semiminor  axis  (Mn)  indicates  a  sense 
counterclockwise  for  the  negative  values. 

of  tidal  current  rotation,  which  is  i 

clockwise  for  the  positive  values  and 

Table  4 

Average  difference  ( H)  between  the  observed  and  computed 

cross- 

and  along-Strait  velocity  components 

e. 

A 

n2 

m2 

Si 

Ei 

Cross-Strait  Velocity  (cm/s)  0.21 

1.35 

mm 

0.46 

Along-Strait  Velocity  (cm/s)  0.20 

HhI 1 

m 

2.22 

Hi  I 

■■ 

0.83 

located  near  the  northern  end  of  the  Strait  as 
shown  in  Figs.  4c  and  d.  Vercelli  (1925)  and 
Defant  (1961),  based  on  the  limited  data,  postu¬ 
lated  that,  in  the  Strait,  there  is  an  amphidromic 
point  of  the  M2  tide,  and  this  amphidrome  is 
located  near  Assab  southwest  of  that  suggested  by 
the  model  results.  The  spatial  distribution  of  the 
N2  amplitude  and  phase  near  the  Hanish  Islands 
resembles  that  of  a  degenerate  amphidromic 
system  (Fig.  4e),  while  the  K2  phase  (Fig.  41) 


decreases  from  less  than  150°  near  Perim  Narrows 
to  30°  just  north  of  the  Hanish  Islands.  The  phase 
of  the  diurnal  constituents  vary  little  in  the  major 
part  of  the  Strait,  and  the  model  results  imply  that 
all  diurnal  constituents  have  anticlockwise  amphi¬ 
dromic  systems  located  north  of  the  Bab  el 
Mandab  Strait  (JC)  and  Pi  at  ~15°N;  0\  and  Q\ 
at  -14°  25'N). 

The  amplitudes  and  inclination  of  the  semimajor 
axis  for  four  principle  constituents,  K\,  Oj,  M2, 


Latitude  Latitude  Latitude 
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Longitude  Longitude 


Fig.  4.  Model  coamplitudes  (in  cm;  solid  line)  and  cophases  (deg,  GMT;  dashed  line)  for  the  (a)  K\,  (b)  O i,  (c)  M2,  (d)  S2,  (e)  N2,  and 
(f)  K2  constituents. 
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and  S2  are  displayed  in  Fig.  5  (inclination  arrows 
are  shown  at  a  subsampled  number  of  the  model 
nodes  for  clarity).  For  all  eight  tidal  components, 
distributions  of  the  semimajor  axis  amplitudes 
show  amplification  of  tidal  currents  as  they  flow 
into  the  Bab  el  Mandab  Strait  from  the  Gulf  of 
Aden.  The  strongest  currents  are  generated  by  the 
K\  (semimajor  axis  amplitudes  of  the  Oh  Pl}  and 
<2i  are,  on  average,  55%,  38%,  and  10%, 


respectively,  of  those  associated  with  the  K\)  and 
M2  (semimajor  axis  amplitudes  of  the  S2,  N2,  and 
K2  are,  on  average,  41%,  25%,  and  5%,  respec¬ 
tively,  of  those  associated  with  the  M2).  Further¬ 
more,  the  amplitudes  of  the  K\  and  M2  currents 
are  comparable.  The  most  energetic  flow  is  present 
in  the  narrowest  part  of  the  Strait  (between  Perim 
Narrows  and  Assab-Mocha  line)  with  the  max¬ 
imum  modeled  amplitudes  found  near  the  Perim 
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Fig.  5.  Model  amplitudes  (cm/s;  a  color  scale)  and  direction  (arrows)  of  the  current  semimajor  axis  for  the  (a)  K\,  (b)  Oy,  (c)  M2,  and 
(d)  S2  constituents. 
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Island  where,  for  instance,  the  speeds  of  the  K}  and 
M2  currents  may  reach  40  cm/s.  Farther  north,  the 
current  amplitudes  associated  with  a  single  con¬ 
stituent  rarely  exceed  20  cm/s.  For  all  constituents, 
the  direction  of  the  maximum  flow  is  generally 
aligned  with  the  along-strait  direction.  Similar  to 
the  results  obtained  from  the  observational  analy¬ 
sis  of  Jarosz  (2002),  the  tidal  currents  are  nearly 
rectilinear  in  the  Strait  (amplitudes  of  the  semi¬ 
minor  axis  are  always  smaller  then  those  of  the 
semimajor  one  and  they  are  not  larger  than  5  cm/ 
s).  The  phase  distributions  (not  shown)  indicate 


that  for  the  same  constituent,  the  currents  turn 
earlier  in  the  shallow  waters  than  those  observed  in 
the  deeper  parts  of  the  Strait.  In  addition,  the 
phase  difference  between  maximum  tidal  velocity 
and  elevation  is  approximately  90°  in  the  central 
part  of  the  Bab  el  Mandab  Strait  (depths  below 
60  m)  for  all  considered  constituents. 

The  next  two  figures  (Figs.  6  and  7)  display,  as 
an  example,  the  modeled  tidal  elevations  and 
currents  observed  during  spring  tides  in  the  Strait. 
Fortnightly  variations  in  this  Strait  are  generated 
not  only  by  interactions  between  the  M2  and  S2 


Fig.  6.  Distribution  of  the  spring  tide  water  elevation  (in  cm)  in  the  Bal  el  Mandab  Strait  at  time  of  (a)  lower  high  water  (LHW),  (b) 
higher  low  water  (HLW),  (c)  higher  high  water  (HHW),  and  (d)  lower  low  water  (LLW)  at  Perim  Narrows. 
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Fig.  7.  Spring  tide  currents  in  the  Bab  el  Mandab  Strait  (a)  ~4  h  after  LHW,  (b)  ~3  h  after  HLW,  (c)  ~3  h  after  HHW,  and  (d)  ~4  h 
after  LLW  at  Perim  Narrows. 


components,  but  also  by  interactions  between  the 
K\  and  0\  constituents;  for  convenience  of  this 
discussion,  we  refer  to  the  largest  tides  generated 
by  all  considered  constituents  in  the  fortnightly 
cycle  as  spring  tides  as  opposed  to  the  lowest  tides 
that  we  refer  to  as  neap  tides.  For  all  displayed 
plots  the  times  of  the  high  or  low  waters  are 
referenced  to  those  observed  at  Perim  Narrows. 

The  spring  tides  in  this  Strait  are  characterized 
by  two  high  and  two  low  waters.  Two  nearly  equal 
high  and  low  water  stages  are  mainly  observed 
near  the  Hanish  Islands  where  the  tidal  elevation 


fluctuations  are  dominated  by  the  semidiurnal 
tides.  South  of  the  islands,  these  maxima  and 
minima  are  more  and  more  dissimilar  due  to  the 
changing  character  of  the  tides  from  semidiurnal 
to  mixed.  The  distributions  of  the  tidal  elevation 
for  both  high  water  stages  of  the  spring  tides  (Figs. 
6a  and  c)  are  fairly  similar.  They  both  show  the 
largest  elevations  in  the  southern  part  of  the  Strait 
with  magnitudes  of  ~70  and  ~90cm  near  Perim 
Narrows  for  lower  high  water  (LHW)  and  higher 
high  water  (HHW),  respectively.  Furthermore,  the 
large  elevation  gradients  are  present  in  the  along- 
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strait  direction  while  simultaneously  showing  little 
or  no  variability  in  the  cross-strait  direction. 
Similar  to  the  high  water  stages,  the  low  water 
elevation  distributions  show  little  variability  in  the 
cross-strait  direction  and  large  variations  in  the 
along-strait  direction;  largest  elevation  gradients 
are  present  during  lower  low  water  (LLW)  and 
much  smaller  ones  occur  at  higher  low  water 
(HLW).  Furthermore,  when  LLW  is  observed  at 
Perim  Narrows  (Fig.  6d),  the  low  water  (negative 
values)  is  generally  present  throughout  the  entire 
Strait,  while  for  the  HLW  stage  (Fig.  6b),  only  the 
part  of  the  Strait  between  Perim  Narrows  and  the 
Assab-Mocha  line  is  actually  in  the  low  water 
stage.  Farther  north,  the  tidal  elevation  distribu¬ 
tion  implies  rather  high  water  stage. 

The  currents  are  in  the  flood  stage  about  4h 
after  LHW  at  Perim  (Fig.  7a).  The  direction  of  the 
flow  shows  some  variability,  but  in  general,  the 
barotropic  currents  flow  into  the  Strait  with  the 
speeds  rarely  exceeding  20  cm/s.  Fig.  7b  displays 
the  speed  and  direction  approximately  3h  after 
HLW  at  Perim  Narrows.  Similar  to  the  previous 
stage,  the  speeds  are  weak  and  rarely  reach  15  cm/ 
s.  The  flow  direction  varies  greatly  and  the  well- 
defined  outflow  (speed  between  10  and  15  cm/s)  is 
generally  confined  to  the  deep  channel.  About  3  h 
after  HHW,  the  tidal  currents  are  at  the  maximum 
flood  stage  (Fig.  7c).  The  flow  is  very  energetic  in 
the  entire  Strait  with  the  maximum  velocity 
reaching  above  100  cm/s  at  Perim  Narrows.  The 
direction  of  the  flood  flow  is  nearly  uniform 
showing  the  well-defined  inflow  across  the  entire 
Strait.  Finally,  about  4h  after  LLW  at  Perim 
Narrows  (Fig.  7d),  the  currents  are  at  the 
maximum  outflow.  The  direction  of  the  flow  varies 
very  little  and  the  outflow  is  observed  across  the 
entire  Strait  with  the  speeds  comparable  to  those 
of  the  maximum  flood  currents. 

During  24  h  of  neap  tides,  one  may  see  one  low 
and  one  high  water  or,  similar  to  the  spring  tides, 
two  high  and  two  low  waters  (data  not  shown). 
The  general  features  of  the  elevation  distribution, 
such  as  higher  magnitudes  near  Perim  Narrows, 
large  gradients  in  the  along-strait  direction,  little 
variability  in  the  cross-shore  direction,  are  iden¬ 
tical  to  those  observed  during  the  high  and  low 
water  of  the  spring  tides.  The  only  difference,  as 


expected,  is  a  magnitude  of  the  tidal  elevations, 
which  is  lower  in  the  Strait  for  the  neap  tides. 
During  these  tides,  the  maximum  speeds  of  the 
flood/ebb  flow  associated  with  HHW/LLW  are 
generally  between  25  and  35  cm/s  south  of  the 
Assab-Mocha  line.  Farther  north,  they  do  not 
exceed  20  cm/s.  Similar  to  the  HHW  flood  or  LLW 
ebb  of  the  spring  tides,  there  is  little  variability  in 
the  direction  of  these  tidal  flows,  which  generally 
show  inflow  and  outflow,  respectively,  across  the 
entire  Strait.  If  the  second  flood  and  ebb  are 
observed,  the  currents  are  even  weaker,  and  their 
direction  is  quite  variable. 


5.  Tidal  residual  circulation  in  the  Strait 

Due  to  variable  bathymetry,  the  tidal  flow  may 
generate  residual  currents  through  nonlinear 
interactions  (Zimmerman,  1978).  These  tidal 
residual  currents  often  play  an  important  role  in 
the  local  mean  circulation,  and  in  the  case  of  the 
Bab  el  Mandab  Strait,  they  may  also  contribute  to 
the  transport  of  the  Red  Sea  waters  to  the  Gulf  of 
Aden/Indian  Ocean  system. 

To  extract  the  time-averaged  residual  motion, 
which  is  induced  by  the  eight  primary  tidal 
constituents  (Ki,  Ou  Pu  Qu  M2,  S2,  N2,  and  K2) 
in  the  Strait,  the  computed  currents  were  averaged 
over  a  29-day  period  (Kowalik  and  Polyakov, 
1998),  and  the  results  of  this  averaging  are  shown 
in  Fig.  8.  In  general,  the  residuals  are  weak  with  an 
average  speed  below  1  cm/s.  In  the  deep  channel, 
to  which  the  Red  Sea  water  outflow  is  confined 
(Murray  and  Johns,  1997),  the  residual  current 
speeds  are  even  lower  and  do  not  exceed  0.5  cm/s. 
At  the  Perim  transect  (see  Fig.  1  for  its  location), 
for  instance,  the  transport  generated  by  the 
residual  flow  in  the  deep  channel  (depths  of  50  m 
or  more)  is  directed  towards  the  Gulf  of  Aden  but 
it  is  very  small  (~3.45  x  10~3Sv,  lSv  =  106m3/s) 
when  compared  with  the  volume  transport  of 
0.37  Sv  associated  with  the  mean  circulation 
(Murray  and  Johns,  1997).  The  strongest  residual 
flow  is  generally  found  near  the  islands  and  in  the 
shallow  parts  of  the  Strait  located  between  Perim 
Narrows  and  Assab-Mocha  line.  Near  Perim  or 
the  Hanish  Islands,  for  instance,  the  speed  of  the 
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Fig.  8.  Time-averaged  residual  tidal  flow  in  the  Bab  el  Mandab  Strait  (a  color  scale — speeds  in  cm/s;  arrows — direction). 


residual  currents  can  be  as  high  as  20  cm/s.  These 
high  speeds  observed  near  the  islands  and  head¬ 
lands  may  be  partly  generated  by  sharp  changes  in 
the  coastline  associated  with  the  model  grid  despite 
its  high  resolution.  The  time-averaged  residual 
flow  also  displays  several  small  clockwise  and 
counterclockwise  eddies  such  as  those  located  west 
of  Mocha  or  south-east  of  Ras  Dehaneba  (see 
circulation  patterns  inside  the  black  rectangles  in 
Fig.  8).  A  lack  of  high  resolution  in  situ  observa¬ 
tions  in  these  areas  prohibits  verification  of  the 
existence  of  these  eddies. 

In  addition  to  the  time-averaged  residual 
currents,  the  time-dependent  residual  flow  was 
also  analyzed.  These  residuals  were  extracted  from 
the  model  computed  time  series  by  removing 
motions  with  periods  less  than  2  days.  The 


amplitudes  of  these  currents  (not  shown)  are 
slightly  lower  than  those  of  the  time-averaged 
flow  with  the  maximum  amplitudes  being  below 
15  cm/s.  Similar  to  the  time-averaged  residual  flow, 
the  most  energetic  currents  are  again  observed 
near  the  Perim  and  Hanish  Islands.  Additionally, 
small  eddies  are  present  in  the  time-dependent 
flow,  and  these  eddies  are  better  developed  during 
the  spring  tides.  These  time-dependent  fluctuations 
are  dominated  by  a  fortnightly  cycle.  Spectral 
analysis  indicates  that  most  of  the  energy  at  low 
frequencies  is  indeed  in  the  fortnightly  frequency 
and  some  in  the  monthly  frequency.  None  of  the 
low  frequency  constituents  is  included  in  the 
forcing.  Thus,  in  this  region,  they  arise  through 
interactions  between  those  constituents,  which  are 
used  to  force  the  model.  The  results  indicate  that 
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Fig.  9.  Amplitudes  (in  cm)  of  the  semimajor  axes  of  (a)  the  fluctuations  generated  by  the  A",  and  Ot  interactions  and  (b)  the  Msr 
constituent. 


among  these  long  period  tidal  constituents,  the  Ms f 
( M2  and  S2  interactions)  and  tide  having  a  13.66- 
day  period  (generated  by  K\  and  Ox  nonlinear 
interactions)  are  dominant  components  in  the 
Strait.  Such  K\  and  0\  interactions  have  been  also 
observed  in  other  ocean  regions  as  described  by 
Brink  (1995)  and  Kowalik  and  Polyakov  (1998). 
Fig.  9  shows  the  spatial  structure  of  the  semimajor 
axes  of  the  Msf  and  13.66-day  components  as 
obtained  from  the  harmonic  analyses  of  a  1 -year¬ 
long  model  time  series.  The  amplitudes  of  this  axis 
vary  between  0.01  and  5.30  cm/s  for  the  Msf  and 
between  0.01  and  5.48  cm/s  for  the  13.66-day 
component.  Additionally,  the  amplitudes  of  the 
Msf  constituent  in  the  Strait  are  generally  smaller 
than  those  associated  with  the  oscillations  gener¬ 
ated  by  the  K\  and  Ox  tidal  components. 


6.  Barotropic  energy  flux  and  dissipation 


The  energy  flux  (Kowalik  and  Proshutinsky, 
1993;  Foreman  et  al.,  1995)  is  given  by 


(5) 


where  h  is  the  water  depth,  £  is  the  tidal  elevation, 
p0  is  Representative  seawater  density  ( =  1026  kg/ 
m3),  U  is  the  current  velocity  vector,  and  g  is  the 
acceleration  due  to  gravity.  This  energy  flux 
includes  the  contribution  of  both  kinetic  and 
potential  energy.  In  the  Bab  el  Mandab  Strait, 
the  flux  of  the  kinetic  energy  is  always  smaller  than 
that  for  the  potential  energy,  on  average,  by  four 
orders  of  magnitude.  For  example,  the  flux  of  the 
K\  kinetic  energy  is,  on  average,  7  x  10~2W/m 
while  the  average  flux  of  the  K i  potential  energy  is 
13.9xl02W/m.  The  same  fluxes  for  M2  are 
5  x  10-2  and  8.5  x  102W/m,  respectively. 

The  average  energy  fluxes  over  a  tidal  period  for 
the  Ki  and  M2  constituents  are  shown  in  Fig.  10. 
For  each  constituent,  the  average  energy  flux  was 
computed  by  averaging  the  flux  time  series 
obtained  from  Eq.  (5)  over  the  constituent  period. 
Additionally,  the  elevations  and  currents,  required 
for  the  flux  time  series  computations,  were 
generated  using  the  model  harmonic  constants 
and  a  time  interval  of  1  min.  It  is  very  apparent 
that  the  average  fluxes  associated  with  two 
dominant  tidal  components  display  different 
behavior  in  the  Strait.  The  Kx  constituent  has 
generally  one  source  of  energy  (Fig.  10a),  which  is 
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Fig.  10.  Mean  tidal  energy  flux  per  unit  length  (in  W/m)  for  the  (a)  K\  and  (b)  M2  constituents. 


the  advective  flux  from  the  Gulf  of  Aden  into  the 
Strait.  The  additional  contribution  of  the  K\ 
energy  from  the  Red  Sea  proper  is  negligible 
because  the  fluxes  associated  with  this  source  are 
at  least  two  orders  of  magnitude  smaller  than 
those  from  the  Gulf.  Similar  to  the  K\,  other 
diurnal  constituents  also  have  one  major  source  of 
energy,  which  is  flux  from  the  Gulf  of  Aden. In 
contrast,  there  are  two  comparable  sources  of 
energy  for  the  M2  constituent  (Fig.  10b):  one  is  the 
flux  directed  from  the  Gulf  of  Aden  and  another  is 
the  flux  from  the  Red  Sea  proper.  The  S2  and  K2 
have  two  fairly  comparable  energy  sources  as  that 
for  the  M2.  For  the  N2,  however,  the  computed 
fluxes  from  the  Red  Sea  are,  on  average,  one  order 
of  magnitude  smaller  than  those  from  the  Gulf. 
Additionally,  Fig.  11  shows  the  average  energy 
fluxes  for  the  spring  and  neap  tides  (both  spring 
and  neap  fluxes  are  25-h  averages).  These  spring 
and  neap  fluxes,  which  give  estimates  for  a 
maximum  and  a  minimum  of  the  barotropic  tidal 
fluxes  in  the  Strait  as  generated  by  eight  principle 
constituents,  clearly  indicate  that  the  major  flux  of 


the  tidal  energy  into  the  Strait  is  again  from  the 
Gulf  of  Aden. 

Table  5  lists  tidal  power  fluxes  normal  to  two 
transects  that  are  located  near  Perim  Narrows  and 
the  Hanish  Sill  (see  Fig.  1  for  their  locations).  All 
estimates  are  the  net  fluxes  through  the  transects 
with  the  flux  into  the  Strait  being  positively  defined 
and  the  flux  out  of  the  Strait  being  negatively 
defined.  In  addition  to  the  constituent  fluxes, 
Table  5  lists  fluxes  of  the  spring  and  neap  tides 
(spring/neap).  The  spring  and  neap  values  are  the 
averages  calculated  over  25  h  periods  beginning  at 
00:00  GMT  on  June  12,  1995  and  ending  at  00:00 
on  June  20,  1995,  respectively,  including  all  eight 
constituents. 

The  values  listed  in  Table  5  indicate  that 
regardless  of  the  constituent  or  spring/neap  tide, 
energy  is  generally  advected  to  the  Strait.  The 
outgoing  fluxes  (see,  for  instance,  diurnal  compo¬ 
nent  and  neap  fluxes)  are  observed  only  along  the 
Hanish  transect  and  their  values  are  just  a  small 
fraction  of  those  transmitted  through  the  Perim 
transect.  In  general,  the  majority  of  the  barotropic 
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Fig.  11.  Mean  tidal  energy  fluxes  (25-h  averages;  in  W/m)  for  (a)  spring  and  (b)  neap  tides. 


Table  5 

Net  power  flux  (in  Gigawatts)  normal  to  the  transects  shown  in  Fig.  1 


Transect 

Ox 

Q\ 

m2 

s2 

n2 

k2 

Spring 

Neap 

Perim  (1) 

mm 

0.004 

0.111 

■a 

1.472 

0.288 

Hanish  (2) 

-0.0003 

0.091 

B9 

-0.003 

tidal  energy  is  advected  to  the  Bab  el  Mandab 
Strait  from  the  Gulf  of  Aden  where  both  diurnal 
and  semidiurnal  tides  are  strong.  Since  the  tides  in 
the  Red  Sea  are  weaker  than  those  of  the  Gulf  of 
Aden  and  semidiurnal  in  nature  it  is  not  un¬ 
expected  that  the  contribution  of  energy  from  this 
Sea  would  be  small,  especially  for  the  diurnal 
components.  The  net  energy  fluxes  along  the 
Hanish  transect  are  indeed  small  with  the  excep¬ 
tion  of  those  associated  with  the  M2  and  K2 
constituents,  which  are  comparable  with  the 
energy  input  from  the  Gulf  as  reflected  at  the 
Perim  Narrows  transect  (Table  5).  In  addition,  the 
model  simulation  with  the  equilibrium  tidal 
potential  excluded  from  the  computations  shows 


that  the  K2  energy  flux  from  the  Red  Sea  to  the 
Strait  is  zero  at  the  Hanish  transect,  while  that  of 
the  M2  flux  is  extremely  small  and  equal  to 
1.9039  x  10-5  GW  at  the  same  location.  These 
findings  imply  that  the  M2  and  K2  energy  fluxes 
across  the  Hanish  transect  are  directly  generated  in 
the  Red  Sea  by  the  equilibrium  tides. 

Due  to  the  small  outflowing  fluxes,  almost  all 
tidal  energy  is  dissipated  within  the  Strait.  Energy 
dissipation  for  all  cases  is  estimated  as  a  net  power 
flux  into  a  polygon  encompassed  by  Perim  and 
Hanish  transects  and  the  coastline  of  the  Strait. 
Note  that  only  one  dissipation  mechanism,  bottom 
friction,  is  included  in  this  particular  implementa¬ 
tion  of  the  ADCIRC  model.  Lateral  and  vertical 
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mixing  and  energy  lost  to  the  generation  of 
internal  tides  are  not  included.  When  the  dissipa¬ 
tion  process  is  viewed  in  terms  of  each  constituent, 
more,  on  average,  semidiurnal  energy  ( M2 — 98%, 
S2 — 98%,  N2 — 85%,  and  K2 — 95%)  than  diurnal 
(A)— 90%,  0,-91%,  Px— 92%,  and  2,-92%)  is 
lost  within  the  Strait  (95%  versus  91%  of  that 
which  enters).  A  high  percentage  of  dissipation  is 
also  observed  during  the  chosen  spring  tide  period 
(98%).  The  dissipation  during  neap  tides  is  slightly 
lower  but  still  large,  equal  to  96%.  Within  the 
polygon  encompassed  by  the  chosen  transects,  the 
dissipation  is  estimated  as  0.226,  0.084,  0.03,  and 
0.0037  GW  for  the  diurnal  constituents  K\,  0,,  P,, 
Q]t,  and  0.202,  0.037,  0.015,  and  0.002  GW  for  the 
semidiurnal  constituents  M2,  S2,  N2,  and  K2.  For 
spring  and  neap  tides,  the  dissipation  is  1.6  and 
0.282  GW,  respectively. 

To  investigate  further  the  dissipation  of  the 
barotropic  tides  in  the  Bab  el  Mandab  Strait,  i.e. 
to  identify  areas  with  high  tidal  energy  losses  and 
find  differences,  if  any,  in  a  dissipation  pattern 
among  different  constituents,  a  rate  of  energy 
dissipation  was  computed.  This  rate  is  defined  as 
(Tsimplis  et  al.,  1995;  Munk,  1997): 

T 

d  =  ^j  CiPo(U2+V2)^2dt,  (6) 

o 

where  Ca  is  the  bottom  drag  coefficient  (  =  0.003), 
p0  is  the  seawater  density  (=  1026  kg/m3),  U,  V 
are  the  velocity  components  at  a  given  grid  point, 
and  T  is  the  tidal  period.  Fig.  12  displays,  as  an 
example,  the  spatial  distribution  of  this  rate  due  to 
the  K ,,  M2,  spring  and  neap  tides.  The  spatial 
distributions  of  the  dissipation  rate  demonstrate 
that  there  are  similarities  in  the  dissipation  pattern 
of  the  diurnal  and  semidiurnal  tides.  Furthermore, 
the  dissipation  processes  are  not  uniform  in  the 
Strait.  The  maximum  rate  of  energy  dissipation, 
for  the  Ki  and  M2  tides  displayed  in  Figs.  12a  and 
b,  is  generally  found  in  the  narrowest  part  of  the 
study  area  located  between  Perim  Narrows  and  the 
Assab-Mocha  line,  i.e.  in  the  region  with  the 
maximum  tidal  velocities.  The  same  pattern  holds 
true  for  the  other  tidal  constituents  considered.  In 
this  region  of  the  maximum  dissipation,  the  rate 


reaches  values  of  or  greater  than  0.01  W/m2  for 
both  Kx  and  M2.  Farther  north,  where  the  tidal 
flow  is  weaker,  the  dissipation  rate  is  generally  on 
the  order  of  10-3  W/m2.  The  dissipation  rate 
estimates  for  other  constituents  are  at  least  one 
order  of  magnitude  smaller  than  those  of  the  K\ 
and  M2.  The  spatial  patterns  of  the  dissipation  rate 
associated  with  the  spring  and  neap  tides  (Figs.  12c 
and  d)  are  very  similar  to  those  of  the  diurnal  and 
semidiurnal  constituents,  and  in  case  of  the  neap 
tides,  the  absolute  values  of  the  dissipation  rates 
are  comparable  to  those  of  the  K\  and  M2.  For  the 
spring  tides,  however,  the  dissipation  rate  esti¬ 
mates  are  generally  one  order  of  magnitude  larger 
than  those  associated  with  the  Kx,  M2,  or  neap 
tides. 


7.  Summary  and  conclusions 

Recent  observations  (Jarosz,  2002)  indicate  that 
tidal  motion  is  a  fundamental,  clearly  non- 
negligible,  aspect  of  the  flow  field  in  the  Bab  el 
Mandab.  The  available  measurements,  however, 
are  spatially  dispersed  and  too  sparse  to  examine 
in  detail  the  tidal  dynamics  in  this  Strait.  Thus, 
barotropic  tides  in  the  Bab  el  Mandab  Strait  are 
numerically  simulated  with  the  high  resolution, 
two-dimensional  form  of  the  finite  element  hydro- 
dynamic  model  ADCIRC  (Luettich  et  al.,  1992; 
Westerink  et  al.,  1994)  to  examine  (1)  variability  of 
the  barotropic  tidal  elevations  and  currents,  (2) 
residual  tidal  circulation  and  its  possible  contribu¬ 
tion  to  the  Red  Sea  water  outflow,  and  (3)  to 
estimate  barotropic  energy  fluxes  and  energy 
dissipation  rates.  To  our  knowledge,  this  model 
is  the  first  numerical  attempt  focused  on  tides  in 
this  Strait. 

Comparisons  with  the  available  observations 
indicate  that  the  model  simulates  well  the  diurnal 
tidal  elevations  when  compared  to  the  observed 
values.  Furthermore,  the  comparison  between 
observed  and  computed  tidal  currents  shows  that 
the  tidal  currents  are  also  well  replicated  in  the 
Strait  for  both  semidiurnal  and  diurnal  constitu¬ 
ents.  The  model  performance,  however,  for  the 
semidiurnal  elevations  is  less  satisfactory  mainly 
due  to  the  errors  between  observed  and  computed 
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Fig.  12.  Tidal  energy  dissipation  rate  (W/m2)  for  the  (a)  K\,  (b)  M2,  (c)  spring,  and  (d)  neap  tides. 


phases  in  the  region  of  amphidromic  systems  for 
these  tides  in  the  Strait.  These  discrepancies 
between  observed  and  computed  phases  are  not 
unexpected  in  such  regions  because  small  errors  in 
the  position  of  the  amphidrome  usually  lead  to 
large  phase  errors. 

The  model  computations  indicate  that  the  tidal 
elevation  range  varies  between  approximately  2  m 
near  Perim  Narrows  to  less  than  lm  near  the 
Hanish  Islands.  The  barotropic  tidal  currents 


generated  by  the  eight  primary  constituents  have 
the  largest  amplitudes  between  Perim  Narrows  and 
Assab-Mocha  line  where  their  combined  speeds 
are  above  1  m/s.  Farther  north,  the  current 
amplitudes  are  reduced  by  approximately  50%  of 
those  observed  in  the  narrow  southern  part  of  the 
Strait.  The  amphidromic  points  for  the  M2,  S2,  N2 
in  the  Strait  and  the  amphidromes  for  the  Qh  Ou 
Pi,  Ki,  and  K2  constituents  just  north  of  the 
Hanish  Islands  together  with  nearly  90°  phase 
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difference  between  the  tidal  elevations  and  cur¬ 
rents  imply  that  the  tides  in  the  Strait  have  some 
characteristics  of  standing  waves. 

Residual  circulation  induced  by  the  barotropic 
tidal  currents  is  rather  weak  with  the  speeds  often 
less  than  1  cm/s  in  the  major  part  of  the  Strait,  and 
it  contributes  little  to  the  overall  transport  of  the 
Red  Sea  waters  to  the  Gulf  of  Aden/Indian  Ocean 
system.  This  residual  circulation  consists  of  mean 
(time-averaged)  and  fluctuating  (deviation  from 
the  mean)  components.  The  fluctuating  residual 
flow  is  dominated  by  fortnightly  oscillations, 
which  are  generated  by  the  interactions  of  the  K\ 
and  0\  constituents  as  well  as  the  interaction 
between  the  M2  and  S2  components. 

The  mean  energy  fluxes  are  small  mainly  due  to 
the  approximately  90°  phase  difference  between 
tidal  velocity  and  elevation,  which  is  observed  in 
the  major  part  of  the  Strait;  however,  these  flux 
estimates  should  be  very  useful  for  comparison 
with  future  calculations  of  the  energy  loss  by 
internal  tide  generation  or/and  turbulent  mixing  in 
the  Strait.  Our  model  results  show  that  the  diurnal 
constituents  have  one  major  source  of  energy,  the 
advective  flux  from  the  Gulf  of  Aden.  In  contrast, 
two  comparable  sources  of  energy  are  identified 
for  the  M2,  and  K2  constituents:  one  is  the 
advective  flux  from  the  Gulf  of  Aden  and  another 
is  the  flux  from  the  Red  Sea  proper.  The  latter 
fluxes  for  both  tidal  components  are  directly 
generated  in  the  Red  Sea  by  the  equilibrium  tides. 
For  the  N2  and  S2,  the  major  flux  comes  from  the 
Gulf  of  Aden.  Furthermore,  very  small  fluxes  from 
the  Strait  to  the  adjacent  basins  indicate  that  the 
Strait  is  a  very  dissipative  environment.  Almost  all 
tidal  energy,  85%  or  more  for  all  considered  cases, 
is  lost  within  the  Strait.  The  distribution  of  the  rate 
of  energy  dissipation  due  to  bottom  friction 
implies  that  the  major  area  of  dissipation  is 
located  between  Perim  Narrows  and  the  Assab- 
Mocha  line. 

In  conclusion,  the  model  results  presented  here 
and  analyses  of  the  available  observations  (Ver- 
celli,  1925,  1927;  Siedler,  1969;  Jarosz,  2002)  have 
significantly  improved  our  understanding  of  the 
tidal  frequency  motion  and  its  importance  for  the 
observed  circulation  in  the  Bab  el  Mandab  Strait; 
however,  they  have  also  brought  to  light  several 


topics  that  need  to  be  addressed  in  the  future.  The 
most  urgent  of  these  is  to  develop  an  under¬ 
standing  of  the  baroclinic  tidal  motion  and  how 
this  motion  interacts  with  water  mass  exchange 
between  the  Red  Sea  and  Gulf  of  Aden/Indian 
Ocean.  Future  work  then  will  focus  on  internal 
tides  and  the  application  of  a  baroclinic  model 
would  be  a  useful  tool  to  help  understand  their 
generation,  propagation,  and  dissipation  in  this 
region. 
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amplitudes  of  the  tidal  elevations  &  the  strongest  currents  are  present  in  the  southern  part  of  the  Strait.  Residual  circulation  induced  by  the  barotropic  tides  is  rather 
weak  in  the  major  part  of  the  Strait,  &  its  contribution  to  the  Red  Sea  water  transport  is  small.  The  model  results  also  show  that  barotropic  energy  fluxes  are  not  very 
large  &  their  direction  depends  on  the  constituent.  All  diurnal  &  one  semidiurnal  (Ay  constituent  have  one  major  source  of  energy,  which  is  the  flux  from  the  Gulf  of 
Aden,  while  there  are  two  sources  of  energy  for  the  M:  S,  &  K,  components:  one  from  the  Gulf  &  another  from  the  Red  Sea.  Very  small  energy  fluxes  from  the 
Strait  to  the  adjacent  basins  indicate  that  almost  all  tidal  energy  is  dissipated  within  the  Strait.  The  distribution  of  the  rate  of  energy  dissipation  due  to  the  bottom 
fliction  Implies  that  the  major  area  of  dissipation  is  located  between  Perim  Narrows  &  the  Assab-Mocha  line. 
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